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ABSTRACT 

We present Keck-I MOSFIRE spectroscopy in the Y and H bands of GDN-8231, a massive, compact, 
star-forming galaxy (SFG) at a redshift z ~ 1.7. Its spectrum reveals both Ha and [Nil] emission 
lines and strong Balmer absorption lines. The Ha and Spitzer MIPS 24 /mi fluxes are both weak, 
thus indicating a low star formation rate of SFR <5—10 M Q y _1 . This, added to a relatively 
young age of ~ 700 Myr measured from the absorption lines, provides the first direct evidence for a 
distant galaxy being caught in the act of rapidly shutting down its star formation. Such quenching 
allows GDN-8231 to become a compact, quiescent galaxy, similar to 3 other galaxies in our sample, by 
z ~ 1.5. Moreover, the color profile of GDN-8231 shows a bluer center, consistent with the predictions 
of recent simulations for an early phase of inside-out quenching. Its line-of-sight velocity dispersion for 
the gas, crfos = 4^7 ± ^ km s_1 > 1S near ly 40% smaller than that of its stars, cr* os = 215 ± 35 km s _1 . 
High-resolution hydro-simulations of galaxies explain such apparently colder gas kinematics of up to a 
factor of ~ 1.5 with rotating disks being viewed at different inclinations and/or centrally concentrated 
star-forming regions. A clear prediction is that their compact, quiescent descendants preserve some 
remnant rotation from their star-forming progenitors. 

Subject headings: galaxies: photometry — galaxies: high-redshift 


1. INTRODUCTION 

The formation and structural evolution of the first qui¬ 
escent galaxies at z > 2 has been the subject of consid¬ 
erable discussion in recent years. Since the first papers 
reporting their remarkably compact nature compared to 
quiescent galaxies of similar stellar mass at low redshift, 
many studies have verified their small sizes and charac- 
terized its growth with cosmic time (IDaddi et all 120051: 
Truiillo e t all 120 071 IBuitraeo et all 120081: Cimatti e t al l 

200a Ivan Dokkum et all 120081 iDamianov et all 120091 : 

Cassata et al. 120111: Szomoru et al. 2012t iNewman et al.1 


2012at iCarollo et al.1 120131 : Ivan der Wei et all I2014D . 
However, the mechanisms responsible for the formation 
of such compact objects are still unclear. 

A crucial step forward in understanding these forma¬ 
tion processes is the recent discovery of a population of 
massive (log(M /M^)> 10), compact star-forming galax¬ 
ies at z > 2 dWuvts et all 1201 Hi IBarro et al.l [20131 
l2014al iPatel et all 120131 : iStefanon et al.l 120 1 31) . These 
galaxies have similar structural properties to the com¬ 
pact quiescent population, i.e., spheroidal morphologies, 
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centrally concentrated luminosity profiles and high Sersic 
indices. The similar properties and number densities sug¬ 
gest that these compact star-forming galaxies (SFGs) 
are the immediate progenit ors of t he s imilarly com¬ 
pact first quiescent galaxies (IBarro et al.l I2013h . Com¬ 
pact SFGs are typically found in a dusty star-forming 
phase characterized by bright far-IR and sub-mm detec¬ 
tions, and relatively normal star formation rates (SFRs) 
similar to those of other star-forming galaxies of the 


as the star-forming main seq 

uence 

(Noeske et al.l 

2007; 

Elbaz et al.||2007t Salim et al|2007i; Pannella et al. 

2009: 

Magdis et al. 2010; Wuvts et al.l 

2011a|: Elbaz et al. 


Pannella et al.l I2014T) . However, they have radically 


different morphologies suggesting that their compact 
nature is the result of strongly dissip ative trans¬ 
formation processes, such as mergers (Hopkins_et_af 
20061 iNaab et al.ll2007t IWuvts et alil2010l:1Wellons et al 


2014) o r accretion-driven dis k i nstabilities (iDekel et al 
200S : iCeverino et al.l 120101 : IDekel &; Burkertl 12014 


Zolotov et al.l I2014T) that funnel a large fraction of their 


gas reservoirs into the center, rapidly building up a dense 
stellar core. 

Additional evidence in support of the evolutionary 
connection between compact SFGs and quiescent galax¬ 
ies came recently when NIR spectroscopy of a sam¬ 
ple of compact SFGs revealed emission l i ne widths of 
a = 200 — 300 km s _1 (|Barro et al.ll2014bl : INelson et al.1 
l20ll . in good agreement with the observed stellar ve¬ 
locity dispersio ns of compact quiescent galaxies of simi¬ 
lar stellar mass (INewman et al.ll2010l:lvan de Sande et al.1 
120131 : iBezanson et al.ll2013fc iBelli et alJl2014blh However. 
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these measurements are based on different dynamical 
tracers, and measured on disjoint populations. Therefore 
the implied evolutionary connection between them is in¬ 
direct. Some caveats to this evolutionary sequence are, 
for example: 1) the broad emission lines may be driven 
by shoc ks and outflows rather than the gravitation al po¬ 
tential (iNewman et alJl2bl2bt iGenzel et al.l 120 bill , and 
2) if compact SFGs do not quench immediately, their cur¬ 
rent dynamics may have little to do with their eventual 
transition into quiescent galaxies. 

A way forward to address these issues is to study the 
properties of both the gas and the stars on the same 
galaxies. These kinematic properties can be used to test 
whether the gas and the stars are probing the same gravi¬ 
tational potential, and the simultaneous measurement of 
emission and absorption lines can be used to estimate 
the age and star-formation history (SFH) of their stellar 
populations. However, these kind of measurements, are 
still observationally challenging, particularly for galaxies 
at z > 1.5. Firstly, absorption line measurements (e.g., 
the Balmer or the Ca HK lines) require NIR spectro¬ 
graphs and long (> 8 hr) integrations on 8 — 10 m class 
telescopes, and probing the stronger emission lines (e.g., 
HaA6563 A, [Nn]A6584 A, [Oiii]A5007 A, H/3A4861 A) 


on the same galaxies often requires additional observa¬ 
tions in other bands. Secondly, the spectra of galaxies 
are typically dominated by either emission or absorption 
lines, and the phases in which both sets of lines are strong 
are usually short and depend on the SFH of the galaxy. 

Interestingly, several papers have identified a sub¬ 
stantial population of recently quenched (i.e., young) 
compact quiescent galaxies at z > 1.5 which are the 
most promising candidates to have both absorption 
lines from an old, underlying population , and emissio n 


li nes from weak, residu al s tar- formation (Newman et al 


201C; Ivan Dokkum fc Brammeil 120101: iToft et al 
2012: lOnodera et all |2012t van de Sande et al.l 2013 : 
Bezanson et all 120131 : iBedreeal et al.l 120131: iBelli et al l 
201 la - iJT Spectroscopic follow-up of these galaxies 


have revealed prominent Balmer absorption lines and, 
occasionally, emission in [On]. The [On]A3726 + 3729 
doublet, however, is typically unresolved at the res¬ 
olution of these observations, and thus it is not a 
reliable kinematic indicator. Additional observations 
of these galaxies for other emission lines in the Ha or 
[O ill] ranges have so far been unsuccessful. 

In this paper, we present Y and H band spectroscopy 
of a massive (log(M/M 0 ) = 10.75), compact SFG at 
z = 1.674 obtained with the MOSFIRE multi-object in¬ 
frared spectrograph (Mc Lean et al.ll2010t iMcLean et al.l 
1201211 on the Keck-I telescope. The spectra of GDN- 
8231 show emission and absorption lines which sug¬ 
gest that the galaxy is quenching rapidly. We model 
these lines to estimate the kinematics of the gas and 
the stars. In addition, we combine the spectra with 
optical medi um-ba nd p hoto metry from the SHARDS 
survey (iPerez-Gonzalez et al.l l2013h and Hubble Space 
Telescope ( HST) WFC3 grism spectroscopy in G102 
(GO13420; PI=Barro) and G141 (GO11600; PI=Weiner) 
to obtain detailed spectral energy distributions (SEDs) 
for the quenching galaxy and 3 older quiescent galaxies, 
observed as a part of the same MOSFIRE survey. From 
the modeling of their SEDs, we estimate their stellar ages 


and SFHs, and analyze the implications for the evolu¬ 
tionary path from compact SFG to compact quiescent 
galaxy. 

Throughout the paper, we adopt a flat cosmology with 
^m=0.3, Da=0.7 and H 0 = 70 km s -1 Mpc -1 and we 
quote magnitudes in the AB system. 


2. DATA 


2.1. Target Selection 


We s elect spectro s copic t argets from the CAN- 
DELS (iGrogin et al.l 120111 : iKoekemoer et al.l 12011ft 
WFC3/F160W (H band) multi-wavelength catalogs in 
GOODS-N (Barro et al. in prep.). This catalog in¬ 
cludes deep, UV-to-NIR ground-based observations in 
several medium and broad bands, HST photometry in 
9 different bands and Spitzer/ Herschel observations in 
the mid-to-far IR. The source extraction and merged 
multi-wavelength spectral energy distributions (SEDs) 
are measured f ollowi ng t he procedu r es de scribed in 
IGalametz et al.l (12013ft and IGuo et al.l ([2013). In addi¬ 
tion, GOODS-N includes complementary HST/ WFC3 
observations in both G102 and G141 grisms, allow¬ 
ing for continuous wavelength coverage from 0.9 < 
A < 1.7 itm with a resol ution better than R= 130 
(e.g., iBrammer et al.l l2012f ). The grism data are re¬ 
duced using the threedhsQ pipeline. The pipeline 
handles the combination and reduction of the dithered 
exposures, and the extraction of the individual spec¬ 
tra. These are background subtracted and corrected for 
contami nation from the over lapping spectra of nearby 
sources (IBrammer et al.l[2012t Momcheva et al. in prep). 
The grism data are joined at shorter wavelengths by 
the 25 optical_ medi um -bands of the SHARDS survey 
(i? ~ 50; iPerez-Gonzalez et al.ll2013 j. Together, these 
datasets provide remarkable spectral resolution on a 
galaxy-by-galaxy basis, that is uniquely suited for SED- 
fitting analysis (see § 3). The stellar masses, rest- 
frame colors, and photometric redslrifts used for tar- 


;et selection are deriv ed via SED-f itting using EAzY 
^Brammer et al.l 120081) and FAST ([Kriek et al. 12009ft . 
assuming iBruzual fc Chariot] (12003ft stellar population 
synthesis models, a Chabrierl (12003ft in itial mass func¬ 
tion (IMF), and the lCalzetti et all ~ (f2000ft dust extinction 
law. The SFRs are computed by adding the unobscured 
and obscured star formation, traced by the UV and IR 
emissio n, follow ing the method described in lWuvts et al.1 
(l2011ai see also lBarro et al.ll2013l l2014aft 
We select galaxies in a transition stage between the 
star-forming main sequence and the quiescent, red se¬ 
quence using —0.75 <log(sSFR/Gyr _1 )< —0.25. This 
corresponds to galaxies ap proximately ~ 2 a belo w the 
main sequence at z ~ 1.75 ( Whitake r et al.ll2012bft . For 
galaxies with log(M/M 0 ) > 10, this threshold roughly 
corresponds to a SFR > 10 M 0 y -1 , which is the 
5cr d etection level with MOSFIRE in ~lh exposures 
(e.g., I Trump et al.l 120131 : iKriek et~aI1 12014ft. We als o 
impose the compactness criterion of lBarro et al.1 (|2013ft . 
M/rl' 5 = 10.3 M 0 kpc -1 ' 5 , to select galaxies with sim¬ 
ilar structural and morphological properties as the qui¬ 
escent population at that redshift. Figure |T| illustrates 
that the sSFR criterion is consistent with the U — V vs. 


15 http://code.google.com/p/threedhst 
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Fig. 1.— Left panel: UVJ diagram for galaxies at 1.3 < z < 1.8 more massive than log(M/M©)> 10 in the GOODS-N field. The colors 
highlight different populations of star-forming galaxies (blue; log(sSFR/Gyr —1 ) > —0.25), quiescent galaxies (red; log(sSFR/Gyr -1 ) < 
—0.75) and quenching, transition, galaxies (green; —0.25 < log(sSFR/Gyr — *) < —0.75) identified according to their UV+IR SFRs. The 
large markers show the quenching, compact SFG (GDN-8231; dark green) and the 3 quenched galaxies (GDN-2617, 17360, 12632; orange to 
red) observed with MOSFIRE. The HST color images ( zJH ) of the 4 galaxies are shown at the bottom. The location of GDN-8231 in the 
UVJ is consistent with the observed spectral and photometric properties, indicating that it is a weakly star-forming galaxy t ransitioning to a 
quiesc ent phase. The quiescent galaxies fall within the selection region for recently quenched galaxies (left of the dashed line: I Whitaker et al.l 
l2012aD . Right panel: Color images (ACS and WFC3), and composite SED of GDN -8231. The grey circles show the (low-resolution) broad¬ 
band photometry, the cyan markers show the SHARDS medium-band data (R~ 50 1 Perez-Gonzalez et alii 20131) and the light and dark blue 
lines show the HST/ WFC3 G102 and G141 grism spectra. The spectral regions probed b y Y and H band MO SFIRE spectra are indicated 
in red. The green and orange lines show the best-fit stellar population templates from [Pacifici et al.l (120121) at a resolution of R= 50 for 
GDN-8231 and the quiescent galaxy GDN-17360. The 3 sub-panels on the right show the zoom-in around the SHARDS, G102 and G141 
data highlighting spectral features in Mguv? [On], HJ and H/3. 


V — J (hereafter UVJ) rest-frame color criterion that has 
been shown to be very successful in identifying' q uiescent 
galax i es, breaking the age/dust degeneracy (Wuvts et al.l 
l2007t iWilliams et al.ll2010l : IWhitaker et al. 1201 lH . The 
spread of the transition sample along the wedge of the 
UVJ quiescent region is primarily driven by differences 
in the dust reddening and in the stellar population ages, 
which suggests a wide diversity of ex tinction levels and 
SFHs among quenching galaxies fe.g.. lWild et al.|[2014T) . 

We select the most promising candidates for spectro¬ 
scopic follow-up by prioritizing bright galaxies with low 
dust reddening to maximize the signal-to-noise (S/N) 
ratio. This additional restriction preferentially selects 
galaxies with bluer UVJ colors, near the so-called post- 
starburst region o f the UVJ diagram (l e ft of the dashed 
line i n Figure [lj IWhitaker et al.l 120111 : iBczanso n et al.l 
l2013fl . This region encompasses colors that are typi¬ 
cal of recently quenched galaxies, with low extinction 
levels and young stellar ages of < 1 Gyr. Galaxies in 
this region also have, on average, lower stellar masses 
(log(M/M 0 ) < 10.8) than older quiescent galaxies wit h 
redder colors (iNewman et al.ll2013l : iBarro et al.ll2014al) . 

The redshift is restricted to the range 1.3 < z < 1.8 
that places the Ha emission line in the H band and sev¬ 
eral other Balrner lines around the 4000 A break in the 

Y band. The observed mask contains 1 candidate that is 
a transitioning , compact SFG (GDN-8231; green circle) 
and 3 recently quenched galaxies (GDN-12632, GDN- 
17360 and GDN-2617; orange to red circles). GDN-8231 
is relatively massive (log(M/M 0 ) = 10.7) and presents 
low levels of star-formation (SFR = 10 M 0 y -1 ) evi¬ 
denced by a weak detection in MIPS 24 //in (/(24) = 
40 //Jy). The low sSFR and mild extinction (Ay = 


0.3 mag) are consistent with its location in the UVJ di¬ 
agram, and suggests that it is on the brink of becoming 
a young quiescent galaxy, similar to the 3 quiescent tar¬ 
gets in our sample. The photometric and spectroscopic 
galaxy properties of the four galaxies are summarized in 
Table 1. 

2.2. Spectroscopic data 

Data were collected on the nights of 2014 April 17 and 
May 11 using the MOSFIRE instrument dMcLean et al.1 
1201011201211 on the Keck-I telescope. The sky conditions 
were clear and the median seeing ranged from 0"5 — 0"7 
seeing. We observed 1 mask configuration in both the 
Y (0.97 < A < 1.12 /rm) and H bands (1.46 < A < 
1.81 /im), with individual exposure times of 180 s and 
120 s, for a total 5.5 h and 2 h, respectively. We used 
2-point dithers separated by 1"5 and slit widths of 0"7. 
The instrumental resolution of MOSFIRE with 0"7 slit 
widths is approximately R = 3200 (~5 A per resolution 
element). The 2D spectra were reduced, sky subtracted, 
wavelength calibrated, and one-dimensionally extracted 
using the public MOSFIRE data reduction pipeline^- 

3. ANALYSIS 

3.1. Kinematic measurements and line ratios 

The 4 galaxies are NIR bright (H < 22 mag) and 
present clear continuum detections (S/N> 5) in both 
the Y and H band spectra. GDN-8231 exhibits multiple 
Balrner absorption lines in the Y band spectrum, from 
Hd to H10, and it is the only galaxy showing Ha and 
[Nil] emission lines in the H band spectrum (Figure 0 . 

11 http://www2.keck.hawaii.edu/inst/mosfire/drp.html 
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Fig. 2.— Left panel: MOSFIRE Y band spectrum of the quenching, compact SFG (GDN-8231) with detected absorption lines. The 
best-fit model used to determine the stellar velocity dispersion and the absorption indices is shown in red. Right panel: MOSFIRE H band 
spectrum of GDN-8231 showing the spectral range around the Her and [Nil] lines. The regions contaminated by sky lines are shown in 
gray. The blue Gaussians show the best fit to the Her (c LOS = 97 km s~ 1 ) and [Nil] (rr LOS = 127 km s~ 1 ) lines using the average continuum 
level. The red Gaussians show the best-fit to the Her line (r^ os = 120 km s —1 ) correcting for the Her stellar absorption as determined from 
the fit to the Y band spectrum. 


The other 3 quiescent galaxies also have clear Balmer 
absorption lines in the Y band but present no emission 
lines in the H band. 

We measure the line-of-sight (LOS) stellar velocity dis¬ 
persion usin g the penalized p ixel-fitting software pPXF 
(|Cappellari fc Emselleml I2004T) . This software fits the 
galaxy spectrum with a model created as a linear com¬ 
bination of the stellar templates that best reproduce the 
galaxy spectrum allowing different weights for each tem¬ 
plates. The stellar templa tes used are t he stars from the 
MILES stellar library (ISanchez-Blazauez et all I2006bt 
ICenarro et al.l 120071) . Before fitting the galaxy spec¬ 
trum, we mask the regions contaminated by telluric at¬ 
mospheric bands and pixels with strong residuals from 
the sky lines, and convolve it with a Gaussian function 
whose width is the quadratic difference between the res¬ 
olution of the MILES stellar library (FWHM = 2.5 A; 
iFalcon-Barroso et al.ll20lH) and the instrumental resolu¬ 
tion of our observations measured in sky lines (FWHM = 
2.7 A). We estimate the uncertainty in the stellar veloc¬ 
ity dispersion by running 1000 Monte Carlo simulations. 
In each simulation, the flux of the spectrum is perturbed 
within a Gaussian function whose width is the uncer¬ 
tainty in the flux obtained in the reduction process. The 
stellar velocity dispersion is measured in each simula¬ 
tion and the uncertainty is the standard deviation of the 
Gaussian distribution of all 1000 velocity dispersion esti¬ 
mates. Based on this method, we obtain a stellar velocity 
dispersion for GDN-8231 of cr* os = 215 ± 35 km s -1 (see 
Table 1 for the quiescent galaxies). 

We measure the LOS gas velocity dispersion, cr®®®, of 
GDN-8231 by fitting a Gaussian profile to the emission 
lines, measuring its FWHM, and subtracting the instru¬ 
mental broadening in quadrature from the FWHM. The 
velocity dispersion is then the corrected FWHM divided 
by 2.355. We fit the Ha and [Nil] lines independently. 
As shown in Figure 0 [Nil] is detected at higher S/N 
ratio because the Ha line is partially contaminated by 
a skyline, and it appears to be self-absorbed (i.e., the 
continuum emission is affected by Balmer absorption). 


As a result, if we fit the lines using the same contin¬ 
uum level the velocity dispersion inferred from Ha is 
smaller than that from [Nil], cr®i)®= 90 ± 18 km s -1 and 
127 ± 32 km s -1 , respectively. Although it is possible 
for Balmer and forbidden lines to have different widths if 
they originate in different regions, the most likely expla¬ 
nation for such a large difference is the self-absorption 
in Ha. To account for that effect, we use the best-fit 
stellar template to the absorption spectra to establish 
the continuum level for the Ha emission, and we recom¬ 
pute the fit. With this method, the inferred velocity 
dispersion is cr®^® = 120 ±34 km s -1 , consistent with the 
[Nil] result. Although the two measurements agree after 
accounting for Balmer absorption, we adopt the higher- 
S/N [Nil] measurement as the more reliable tracer of 
gas velocity dispersion. 

The spectra are not flux calibrated. However, we 
can use the equivalent width of Ha corrected for stel¬ 
lar absorption (EW(Ha) corr = 8.3 ± 0.7) and the con¬ 
tinuum flux, inferred from the SED modeling, to cal¬ 
culate the H a line flu x and S FR. Using the empirical 
relation from iKennicutt ( 199 8) and the attenuation in¬ 
ferred from SED-fitting (Ay = 0.3), we obtain values 
of SFR = 3 — 6 M 0 y -1 , depending on the extra neb¬ 
ular extinction with re spect to the continuum (Af^ = 
2.44 — 1.86 A cont ; e.g., iCalzetti et al.l[200(1 iPrice et all 
[2on . These values are consistent with the estimate from 
MIPS 24 fim data. The Ha line flux corrected for stellar 
absorption also allows us to estimate the intrinsic value 
of the line ratio [Nn]/Ha= 1.2. Even without measur¬ 
ing the [O ill]/H /3 line ratio to full y constrai n the ioniza¬ 
tion diagnostic diagram (i.e., BPT: lBaldwin et aT1ll98ll) . 
an [Nil]/Ha ratio of the order of unity already suggests 
that the nebular emission is at least partially powered 
by an AGN. The ga laxy, however, is not detected in 
the 2Ms X-ray data ([Alexander et all12003 ) which im¬ 
plies that the AGN is relatively weak, perhaps shut¬ 
ting down along with the star formation in the galaxy. 
This result, together with the large fraction of AGNs 
(~ 40%) found in compact SFGs at z = 2—3 (jBarro et al.l 
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l2014af) . suggest that the black hole growth at high red- 
shift closely follows the star-formation history (i.e., it 
starts and shuts down with the star-formation), as op¬ 
posed to the observed trend in the local Universe, where 
there seems be to a delay between the starburst and the 
peak of the AGN phase dWild et, al.l 1201(1 : lYesuf et, ahl 
12014 iHernan- Caballero et al.l I2014T) . In that case, AGN 
may play a more active role in the quenching of star for¬ 
mation at z ~ 2. 


3.2. Luminosity profile and Morphology 

The MOSFIRE spectra of GDN-8231 are spatially un¬ 
resolved and thus do not provide additional insights 
on the kinematic profile of the gas and the stars, or 
the spatial distribution of star formation. The high- 
resolution HST imaging, however can be used to de¬ 
termine the overall structural properties and the radial 
distribution of star-formation in the galaxy. Figure [3] 
shows the surface brightness and color profiles of GDN- 
8231 computed from the fitting of the HST-based SEDs 
measured at each radius. The profile shows a positive 
color gradient of d(NUV — V)/dr = 6 • 10” 3 mag/kpc 
that is indicated by a bluer center with a 20% smaller 
r e in the rest-frame NUV with respect to the V band 
(r et NUV/fe,v = 0.8). Assuming that the NUV is a good 
tracer of the SFR, these measurements suggest that the 
SFR profile is more centrally-concentrated than the stel¬ 
lar profile. A possible interpretation for this result is that 
we are witnessing the quenching of a central starburst 
right after it has formed a dense stellar core (Ei k pc = 
9.7 M 0 kpc” 2 ), which seems to be a pre-requisite for 
quiescent galaxies (Cheung et al.ll2012l iFang et all 1 20 1.4 


Ivan Dokkum et all 20141 . Such interpretation is an ex¬ 

cellent match for the predictions of simulations that de¬ 
scribe the formation of compact galaxies as the result 
of strongly dissipative gas-rich events, such as merg¬ 
ers and/or disk inst a bilities (e.g., Ibekel fc Burkertll20l4 
iZolotov et al.l 12014 iWellons et all 120141 . These pro¬ 
cesses are characterized by a wet-inflow (i.e., gas-inflow 
rate > SFR) that builds-up the central gas density, 
thereby enhancing the SFR at the center and growing 
a dense stellar core. The weakening of such inflow marks 
the onset of inside-out quenching due to gas depletion. 
In that onset phase, the SFR is still high and thus the 
color gradient is still positive, as observed in GDN-8231. 
However, as inside-out quenching progresses, the SFR 
profile flattens, and the color gradient turns negat iv e, as 
o bserved in compact qu iescent galaxies (|Guo et, al.ll2012l : 
?: lSzomoru et al1l2012H . 

GDN-8231 has high-Sersic values consistent with 
ot her com p act st ar-forming and quiescent galaxies 
(iBarro et all 12013 1. Its visual appearance in the H 
band is smooth and spheroidal. However, the ACS im¬ 
ages, that probe the rest-frame UV, show small irregular 
patches, perhaps indicative of its prior, diss ipative c om¬ 
pactio n event. Adopting the assumptions of lMiller et al.1 
(1201 11 1, we infer a viewing-angle inclination of i = 42° 
from an axis-ratio value of 6/a = 0.74. 


3.3. Spectral indices and SED fitting 

One of the main obstacles to estimating stellar 
population properties from the analysis of SEDs is 
that at the typical resolution of broad-band surveys 



radius [kpc] 


Fig. 3.— Rest-frame surface brightness profile and color profile of 
GDN-8231. The rest-frame profiles are computed by interpolating 
at each radius the best-fit SED derived from the observed surface 
bri ghtness profi le in 9 HST bands measured with IRAF/ellipse (see 
ILiu et ahl 120131 and Liu et al. in prep for more details). At z = 
1.67, the rest-frame NUV and V bands roughly correspond with 
the observed i and H bands, respectively. The arrows indicate 
the effective radius in those bands from the best-fit Sersic profiles 
obtained using GALFIT. The grey line indicates the PSF half width 
at half-maximum (HWHM) in the H band. The insets show the 
images of GDN-8231 in the i (PSF-matched to H band) and H 
bands. The black circle has radius of 1"(~ 8.4 kpc). GDN-8231 
has an positive color gradient and relatively low sSFR, which are 
consistent with the expectations for an early phase of inside-out 
quenching after a central starburst. 


(FWHM~0.2 /rm/(l-|-z); R~6) the most relevant con¬ 
tinuum and emission/absorption line features are usu¬ 
ally diluted, whic h results in large uncertainties in in¬ 
ferred properties (jMuzzin et al.1 120091 : iGonrov fc Gunnl 
IM3). One way to circumvent this problem is by 
using higher spectral resolution data to obtain more 
accurate measurements of line strengths and spectral 
indices, which are key indicators of stella r age, and 
present /past star-formation activity (e.g., Kelson et al.1 
120011: iKauffmann et al.ll2003lf . 

From the MOSFIR E Y band spectra we measure 
the H6 a Lick index (IWorthev fo Ottavianilll997 1 using 
PPXF to estimate the best-fit value and uncertainty. 
This index h as been trad itionally used as a n age an d SF H 
indica tor dTrager et al.l 120001 : iSanchez-Blazauez et al.l 
I2006al) . From the H band spectra we measure the 
Ha Equivalent Width, EW(Ha), which is sensitive to 
the instantaneous (last few Myr) SFR. Only GDN-8231 
presents Ha emission, the other 3 quiescent galaxies 
have strong continuum detections that place reliable up¬ 
per limits on EW(Ha). In addition to the MOSFIRE 
spectra, we use the SHARDS medium-band data and 
the HST grism data to probe for additional absorp¬ 
tion features and continuum breaks. The right panel 
of Figure |T] illustrates the high spectral resolution (~ 
10 — 100 x better than broad-band filters) of the merged 
medium-band/grism SED, revealing Mgn absorption at 


2800 A in SHARDS and Balmer absorption lines and 
the 4000 A break in G102 and G141. We quantify the 
later using the D4000 index (e.g., iBalogh et al.1 119991 : 
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Age of the Universe [Gyr] 



Fig. 4.— Evolutionary tracks in sSFR (left) and luminosity-weighted age (right) vs. age of the Universe for different SFHs. The 
blue shaded region depicts the star-forming main sequence determined from the average SF H of SF Gs drawn from the model library of 
IPacifici et al.l (|2012) . This region agrees well with the observational results of I Whitaker et all (|2014fl for SFGs of intermediate mass. The 
black lines illustrate the evolution of 2 galaxies that have a secular growth from z ~ 4 to 2 = 1.7, followed by either fast (solid) or slow 
(dashed) quenching of star formation. In a fast quenching galaxy the luminosity-weighed age grows linearly with time (passive evolution). 
However, in a slow quenching (or main sequence) galaxy the luminosity-weighed age increases more slowly (i.e., the slope is < 1). At 
model can describe either a fast or slow quenching. However, a short r (gray line) would provide unrealistic results for galaxies described 
by a two-phase SFH (e.g., a main sequence + fast quenching). 


IKauffmann et al1l2003f) . The G102 grism shows also a 
weak [Oil] emission line (EW 15 ± 5 A). However, there 
is no clear sign of [O ill] emission in G141. The later is 
expected to have also low EW < 10 A, and thus can 
be partially hidden by the H/3 absorption, as hinted 
in the stacked s pectra of recently quenched galaxies in 
iWhitaker et al.1 (120131) . 

The spectral indices are usually analyzed by compar¬ 
ing measurements to a grid of values derived from stellar 
population synthesis models. Here we follow a slightly 
different approach to combine information from both 
the indices and the overall UV-to-N IR SED bv using 
of the SED-fitting code of IPacifici et al.1 (12012L here¬ 
after P12). The code performs a simultaneous fitting of 
the low and high resolution data and includes priors on 
the EW of emission and absorption lines to obtain bet¬ 
ter constraints on the SFR ar id S FH of the galaxy ( see 
also IPacifici et al.ll20l4120151 and iBarro et al’ll2014a] for 
more details). The galaxy templates are computed from 
non-parametric SFHs adapted from semi-analytic mod¬ 
els (SAMs) and include both the stellar continuum and 
the nebular emission. The SAM SFHs provide a richer 
parameter space featuring short-timescale variations of 
the SFR (burst and truncations) that are missing in ex¬ 
ponentially declini ng ( exp[-t/r]) or delayed (txexp[-t/r]) 
models (see also § I4TI) . 

4. RESULTS 

4.1. Stellar ages and SFHs 

The stellar age and SFH are not independent proper¬ 
ties and thus, the chosen paranretrization of the latter 
(single burst, N-bursts, r models, SAMs, etc) often leads 
to strong degeneracies in the age and the characteris¬ 
tic timescale(s) of star formation (e.g., age-r; iConrovI 


120131 ). These degeneracies, however, can be reduced us¬ 
ing higher spectral-resolution datasets to probe features 
that are sensitiv e to different star -formation and quench¬ 
ing t i mescales feriek et al.l 1201 ll: iPerez-G onzalc z et al.l 
120131: IBelli et a,l .11201 51) . 

Figure [4] illustrates the implications of assuming dif¬ 
ferent SFHs to estimate the age of the galaxies. The 
blue shaded region shows the evolution of the sSFR and 
luminosity-weighted age vs. time for galaxies in the star¬ 
forming main seque nce, as describ ed by the SAM SFHs in 
the model library of lPacifici et al. ( 201 2). Those galaxies 
are thoug ht to be growing in a relatively smooth, se cu- 
lar mode (lElbaz et, al.l 120071 : iRodighiero et al.1 120101 ) in 
which gas in flow and S FR have reached a steady-state 
phase fe.g.. iDekel et al.ll2013l) . The main sequence pic¬ 
tured by S AM SFHs follows clos ely the observational re¬ 
sults fe.g.. IWhitaker et al.ll2014D which are also in good 
agre ement with the predictions of other em pirical models 
(e.g., iPeng et aLll2010t iBehroozi et al.1120131 ). In the main 
sequence paradigm, quenching can be interpreted as the 
departure from a stable growth phase, which can be ei¬ 
ther fast (e.g., SFR truncation; solid line) or slow (shal¬ 
lower slope; dashed line). After fast quenching, age in¬ 
creases linearly with time (i.e., passive evolution), while 
for slower quenching the slope is shallower, and the age 
increases slowly due to slowly declining star-formation 
(right panel of Figure [4]). In r models, the timescale con¬ 
trols the quenching time, and it can be tuned to describe 
fast or slow quenching (grey line in Figure [4]). However, 
as one parameter characterizes the whole SFH, the pre¬ 
diction is not realistic for galaxies that are described by 2 
or more distinct phases, such as a secular growth followed 
by fast quenching. 

Panels a) and b) of Figure[5]show index-index diagrams 
sensitive to the SFH and quenching time. The colored 
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redshift 



Age of the Universe (Gyr) 
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Fig. 5.— Panels a) & b): Index-index diagrams for the galaxies in the sample. The circles show the values measured in the spectra. The 
colors are the same as in Figure[l] The colored lines depict the evolutionary tracks for the best-fit SAM SFHs. The tracks start at the onset 
of quenching (colored squares) and continue for ~1 Gyr after the redshift of observation. The blue shaded region and the dashed line show 
the evolution of the star-forming main sequence and the slow quenching galaxy depicted in Figure |4] and Panels c) and d). The Ha and 
H5 indices are sensitive to recent changes in the SFR, while D4000 is a good tracer of the SFH averaged over longer timescales. In a fast 
quenching, the EW(Ha) drops abruptly while H5 a increases due to the appearance of prominent Balmer absorption lines. In contrast, a 
slow quenching galaxy has higher EW(Ha) and weaker H^a for the same value of the D4000 index. Panels c) & d): Same as Figure [4] for 
the galaxies in our sample. The grey stars show the age of the oldest and youngest quiescent galaxies at z ~ 2 in I Whitaker et aEll2013l All 
4 galaxies are consistent with fast quenching. However, the details of the quenching path lead to different luminosity-weighted ages. For 
example, GDN-12632 was the first of quiescent galaxies to quench and is the youngest, while GDN-2617 was the last to quench and it is 
the oldest. 


circles indicate the values measured in the MOSFIRE 
spectra of the 4 galaxies, and the lines show the evolu¬ 
tionary tracks of their best-fit SAM SFHs. The D4000 
index is a good tracer of the average age of the stellar 
populations, while the EW(Ha) and Hi5a are more sen¬ 
sitive to recent star-formation. For slow quenching, the 
EW(Ha) in emission decreases slowly while the Balmer 
absorption is relatively small and constant with time, 
H(5 a ~ 1.8. In contrast, fast quenching causes a rapid 
decline in the Ha emission followed by an absorption 
plateau with EW(Ha)~ —2 A, and a strong H5 a absorp¬ 
tion. This relatively short-lived phase (~1 Gyr), charac¬ 
terized by the lack of Ha emission and strong absorption 
in other Balmer lines, typical of A-stars, is know n as a 
post-starburst phase (pSB; e.g.. lWild et alJ 1201 0( 1. 


The spectral properties of GDN-8231 agree well with 
the sSFR selection criteria indicating that the galaxy 
is on a fast quenching path. Besides the Ha emis¬ 
sion, which indicates the presence of weak, ongoing star- 
formation, the SED fit, the D4000 index and the H(5 a val¬ 
ues suggest that the emission is rapidly declining. In fact, 
given the short quenching time, the chances of finding a 
galaxy in this phase are small, which makes GDN-8231 
quite unique. The 3 quiescent galaxies in our sample 
are further along in their evolution, as indicated by their 
stronger D4000 ~ 1.6. However their indices and best- 
fit SFHs are also consistent with fast quenching. The 
similar SFHs are another indication that GDN-8231 is 
evolutionarily linked to the quenched galaxies. For a def¬ 
inition of quenching time, t q , as the elapsed time from 





















having EW(Ha) values consistent with those of a main 
sequence galaxy, to EW(Ha)~ 0 (tracks in panel-a), the 
values range from a few Myrs for GDN-8231 (i.e., SFR 
truncation) to 500 Myr for GDN-17360. Interestingly the 
SAM SFHs suggest that galaxies show pSB features as a 
result of fast quenching from the star-forming main se¬ 
quence, without being preceded by a strong starburst 
(i.e., ~ 3x higher SFR than the main sequence). A 
more practical definition of t q as the elapsed time from 
log(sSFR/Gyr _1 )= 0 to -1 (tracks in panel c)), results 
in similar values of t q = 300 — 800 Myr. 

The fast quenching scenario for all 4 galaxies 
agrees well with recent works that find strong B aimer 
absorptions on similarly young qu ie scent galaxies 
at z ~ 1.5 dNewman et all l201Ct lOnodera et al.l 


201^__ |Bezanson_et all 1 20131 Ivan de Sande et all 1 20131 
BeilLet_aL| 120. 4a ). I t differs, however, with the results 


of lKriek et, al.l (120111 ). who studied a sample of stacked, 
high-resolution SEDs of galaxies at 0.5 < z < 2.0 find¬ 
ing EW(Ha) values consistent with gradually declin¬ 
ing SFHs. Assuming the difference is not caused by 
the slightly different redshift range, a plausible expla¬ 
nation for this discrepancy is that there are in deed dif¬ 
feren t evolutionary quenchin g tracks (e.g.. iMartin et. all 
2007t iGoncalves et al.1 120121 : iSchawinski et al.l 120141 or 


Belli et al.ll2015D . and current spectroscopic samples at 


z > 1.5 are still not fully representative of the whole 
quenching population. 

Panel d) of Figure [5] shows that the luminosity- 
weighted ages of the sample ranges from t w = 700 Myr 
for GDN-8231 to t^, = 1.1 — 1.3 Gyr for the 3 quies¬ 
cent galaxies. The SAM SFHs also suggest a formation 
redshift (i.e., onset of star formation) of ^form > 6, a 
half-mass assembly by z ~ 3 and the onset of quench¬ 
ing by z ~ 2. Interestingly, the latter is similar to the 
Zform of the galaxies inferred from r models, which sug¬ 
gests a rapid build up (r < 100 Myr). This, however, 
is most likely a limitation of single-parameter r mod¬ 
els, which are biased towards short values of r in or¬ 
der to reproduce the strong pSB features in the SED. 
Secular SFHs with rise and decay timescales of several 
hundred Myr appear to be more realistic, a nd are pre¬ 
dicted by physically motivated models (e.g.. iPeng et al.l 
120101: iBehroozi et al.l 120131 : iGladders et all 120131b How- 
ever, obtaining a direct, reliable measurement of the for¬ 
mation timescale would require additional measurements 
such as the metallicity [Z/H] or the a element abun¬ 
dance la/Fel l _wMch_trace_the_evolution of the metals 
(e.g., lConrovir2013l : lOnodera et al.ll20l4l) . 

The stars in panel d) of Figure [5] show the age of the 
oldest, t w = 1.5 Gyr, and young est, t,„= 0.9 Gvr, qui¬ 
escent g alaxies at z ~ 2 fro m IWhitaker et al.l (|2013l ) 
(see also iNewman et, al.l 120131 ) . The evolutionary tracks 
of the quiescent galaxies in our sample are consistent 
wi th being direct descen dants of the youngest galaxies 
of iWhitaker et afl (120131) . and GDN-8231 reaches a sim¬ 
ilar age by 2 = 1.5. This implies that as new quenching 
galaxies join the red sequence the age spread in the qui¬ 
escent population will increase by roughly 1 Gyr from 
z = 2 to z = 1.5. Interestingly, if all the galaxies in 
our sample follow a 'pure passive evolution since z = 1.7, 
their ages by z = 0.6 would be t™ ~ 4. 5 — 5 Gyr, which is 
roughly 2 Gyr older than the results of lChoi et al.l (|2014D 



Fig. 6.— <t los vs. M* for different galaxy samples. The point up 
and down triangles show emission (gas) and absorption (stars) line 
measure ments, respec t ively. T he blue triangles depic t 13 compact 
SFGs in |5arro et ah (|2014aT) . and one galaxy from [Nelson et al. 
(2014 ). The red triangles sho w a compilation of quie s cent galax¬ 
ies (|van de Sande et al.l 120131 ; IBezanson et alj 120131 ; IBelli et al.l 
t2014alfB) ' at z > 1.5. The overlapping distributions for compact 
SFGs and quiescent galaxies suggests that both populations have 
similar kinematic properties. This is supported by the agreement 
in t he gas and st ellar dispersions of COS-10289 (IBelli et a,I.II2014bl . 
and IBarro et alll2014bl for the stellar and gas kinematics). How¬ 
ever, GDN-8231 (green) has ~ 40% lower dispersion in the gas 
than in the stars. A possible explanation is that the gas has colder 
kinematics than the stars (i^/oy > 1), and thus its line-of-sight 
dispersion is prone to stronger projection effects. The gas and 
stellar dispersion for the simulated galaxy V12 (see also Figure [7} 
illustrates the bias towards lower values of with respect to 

°los f° r v ^ w 'ng-angles closer to face-on (magenta). 

for the oldest quiescent galaxies at that redshift. This 
could be an indication that some quiescent galaxies re¬ 
tain low levels of star formation, or perhaps experience 
minor wet mergers that rejuvenate star formation. 

In summary, the 4 galaxies in our sample have rel¬ 
atively young ages of t u ,< 1 Gyr and present spectral 
features consistent with fast quenching of the star for¬ 
mation. In particular, GDN-8231 seems to be caught in 
an early stage of a rapid truncation. 


4.2. Kinematic properties 

Figure [6] shows the M*—cr LOS relation for a com¬ 
pilation of compact SFGs from iBarro et al.l (l2014aD 
and quiescent galaxies at z > 1.5 from the litera¬ 
ture (Iva n de Sa nde et al l 12013: IB ezanson et, al.l 120131 : 
IBelli et al.l l2014aHbh . In IBarro et all (|2014aF ~ we found 
that 1) the gas kinematics of massive, compact SFGs are 
consistent with the stellar kinematics of compact quies¬ 
cent galaxies, and 2) the dynamical masses of both pop¬ 
ulations are in good agreement with their stellar masses. 
We interpreted those similarities as evidence for an evolu¬ 
tionary connection. In particular, since both populations 
have similar structural, kinematic and dynamical prop¬ 
erties, compact SFGs simply turn into compact quiescent 
galaxies by quenching their star formation. In addition, 
one galaxy (COS-10289), seemed to support this scenario 
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Fig. 7. — Stellar (left) and gas (right) kinematic profiles for the simulated galaxy V12 at z = 2.3 llCeverino et all 120141 : IZolotov et al.l 
[20140. The kinematic properties are measured in cylindrical beams with a depth of 8 kpc. The solid and dashed black lines depict the 
intrinsic rotation and dispersion. The magenta and purple lines show the c LOS profiles for a line-of-sight inclination of i = 90° (edge-on) 
and i = 45°. The triangles show the integrated, mass-weighted values at the r = r e . The line-of-sight dispersion can be written as 
°LOS = /3( s i n * i.y/j) 2 + Cr, where /? depends on the inclination and density profile of the galaxy. This implies that for a rotation dominated 
component the value of <r LOS depends more strongly on projection effects. In V12 the gas has higher rotation than the stars (u^/oy ~ 1.4 
vs. v^/crr ~ 5) and therefore, the ratio of velocity dispersions can be as high as j ~ 1-5 for a LOS inclination of i = 45°. The 

insets in the upper-right show the 10x10 kpc face-on density maps for the stars and the gas. The black circle has r = r e . 


by showing agreement between its gas and stellar kine¬ 
matics. However, the low S/N ratio of the emission line 
measurement and the presence of an X-ray AGN within 
the galaxy place caveats on the interpretation of the kine¬ 
matic properties of the gas. 

Quite surprisingly, the better-constrained measure¬ 
ments for GDN-8231 (green triangles in Figure |6j show 
that the gas has a lower velocity dispersion than the stars 
by a factor of a * Q = 1.7 ± 0.5. Naively, we ex¬ 

pected a value ~ 1 or, given that GDN-8231 is quench¬ 
ing, larger values of cr 8 ^ caused by strong feedback pro¬ 
cesses Jo L g ;j _shocks o^outflows; iDiamond-Stanic et al.1 
1201 21 : iGenzel et alJl2014l ). In turn, the smaller values 
of suggest that the gas is in dynamical equilibrium 
and that it has colder kinematic properties than the stars 
(i.e., higher ty^/ay in the gas). If that is indeed the case, 
the lower dispersion in the gas could be the result of: 
1) widespread star-formation activity in a disk observed 
at low inclination (i.e., close to edge-on), or 2) a cen¬ 
trally concentrated star-forming region probing, on av¬ 
erage, lower values of the rotational velocity v,),, which 
grows inside-out. As described in § 13.21 GDN-8231, has 
an inclination of i = 42°, and 20% smaller effective ra¬ 
dius in the rest-franre NUV, which can lead to the smaller 
values of oi s ^® compared to a* . 

b D b # b ^ ^ # 

In order to provide a better intuition, and quantify how 
much projection effects and/or a concentrated SFR pro¬ 
file affect the measurement of cr LOS , we study the kine¬ 
matic profiles of the gas and the stars in V12, a high- 
resolution, hydrodynamic gal axy sim ulation (~2 5 pc 
grid) drawn from the s ample of IZolotov et al.l (120141) and 
ICeverino et al.l (120141) . As described in IZolotov et all 
1 20141) . these galaxies have similar stellar and structural 
evolution as the compact SFGs and therefore, V12 pro¬ 
vides a excellent proxy for the analysis of the kinematic 
properties of GDN-8231. Figure [7] shows that the stars 
in V12 have comparable rotation and dispersion whereas 


the gas is rotation dominated (stellar iyi/oy~ 1.4 vs. 
gas u</>/oy~ 5). As a result, is more sensitive to 
projection effects and shows smaller values for low in¬ 
clinations. For example, the integrated, mass-weighted 
er 8 at r = r e show ratios of a* . /a 8 ?® ~ 1.0 and 1.5 
for an edge-on (i = 90°) and an intermediate inclination 
(i = 45°), respectively (see values in Figure [G]). There¬ 
fore, we conclude that a large fraction of the observed ra¬ 
tio of velocity dispersions in GDN-8231 can be accounted 
for by projection effects. 

On the other hand, given the relatively flat c 8 /)® pro¬ 
file of V12 (right panel of Figure 0, the change in the 
integrated value of 0^(r e ) would be small for a gas den¬ 
sity (ocSFR) profile more centrally concentrated than the 
stellar profile (i.e., r e ,* > r e ,gas)- Using the ratio of effec¬ 
tive radii in GDN-8231 (r e ^uv/ r e,v = 0.8) as an exam¬ 
ple of different gas-to-stellar concentrations, we find that 
integrating a 8 ^ only up r = 0.8r e decreases its value 
by less than 5% for the edge-on case. Nevertheless, the 
effect of the gas density profile in cr 8 ^ can be larger in 
GDN-8231 if: 1) it had a slowly-increasing rotation curve 
and a flat ay in the center, and/or 2) the emission line 
region, traced by Ha, had an even smaller r e than the 
NUV luminosity profile. The latter is more plausible 
given the quenching nature of GDN-8231 and the shorter 
star-formation timescales probed by Ha. 

Lastly, note that Figure [7] shows only the radial com¬ 
ponent of the intrinsic dispersion (ay) for V12. In the 
case of anisotropic dispersion, for example due to strong 
collimated winds in the center of the galaxy (oy oy), 
the observed value of o' 8 /)® for a face-on inclination would 
be much larger than that of o-* os . 

4.2.1. Dynamical mass 

We estimate the dynamical mass of GDN-8231 from 
a * Q which, as described above, is less sensitive to projec- 
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tion effects. Following the virial equation: 

M dyn (r < r e ) = K a ^° ( f e ( 1 ) 

where K depends on the galaxy’s mass distribution, 
the inclination, and velocity field. We use K = 2.5, 
which is the most widely adopted value for stars i n dis¬ 
persion dominated galaxies (e.g., | Newma n et al.1 120101 : 
Ivan de Sande et aj]l2013t IBelli et al.ll2014bf l. and is valid 
under a variety of galaxy geometries and mass dis¬ 
tributions (e.g., iBinnev fo Tremainel 12008 1. We ap¬ 
ply a small correction to scale the observed a* to 
the value i n a circular aperture o f radius r e , oy = 

I. 05xf7* og dCappellari et al.l 120061 : van de Sande et all 
l201ll l. The inferred dynamical mass log(Md yn /MQ) = 

II. 1 is slightly larger than the stellar mass, but still 
consistent within the uncertainties in both the stellar 
mass and velocity dispersion. Previous works have re¬ 
ported a similar offset towards larger dynamical masses 
of log(M d yn/M*) = 0.1 — 0.16 for recently quenched qui¬ 
escent galaxies z > 1.5 (IBelli et a,Ul20im 

5. CONCLUSIONS 

We present Keck-I MOSFIRE NIR spectroscopy of 
GDN-8231, a massive, compact SFG galaxy at z ~ 1.7. 
This galaxy was selected with sSFR and rest-frame col¬ 
ors matching an intermediate stage between star-forming 
and quiescent. The Y and H band spectra reveal strong 
Balmer absorption lines and Ha and [Nil] in emission. 
The emission and absorption lines yield spectral indices 
and the kinematics of the gas and the stars. 

The spectral indices, SED-modeling, and the compar¬ 
ison to 3 compact quiescent galaxies at similar redshift, 
indicate that GDN-8231 was caught in a rare, early 
stage of fast quenching. Still relatively young, with a 
luminosity-weighted age of 700 ± 250 Myr, GDN-8231 
will mature to become a compact quiescent galaxy by 
redshift z = 1.5. The rapid truncation of the SFR is 
evidenced by the low EW(Ha) and weak MIPS 24/urn 
flux. The color profile is bluer in the center, which is 
consistent with the predictions of recent simulations for 
an early stage of inside-out quenching. The line ratio of 
[Nil]/Ha ~ 1 suggests the presence of a weak (not X-ray 
detected) AGN, a common finding among most compact 
SFGs (Barro et al. 2014a). Using SFHs based on SAMs, 
we find that the assembly of GDN-8231 is consistent with 
an early (z ~ 6) onset of star-formation, a secular build¬ 
up in the star-forming main sequence, forming 50% of its 
stellar mass before z = 3, and fast quenching at z ~ 1.7. 

In Barro et al. (2014b), we found that compact SFGs 
and quiescent galaxies have similar line-of-sight veloc¬ 
ity dispersions for the gas and the stars, suggesting 
similar kinematics in progenitors and descendants. In 
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GDN-8231, however, the dispersion of the gas, cr®= 
147 ± 32 km s -1 , is 40% smaller than that of the stars, 
°los = 215 ± 35 km s -1 . This difference can be explained 
if the gas has colder kinematics (rotation dominated) 
than the stars, and therefore: a) a is smaller if the 
viewing-angle is low (close to face-on), and b) a ^ is 
smaller if the emission-line (star-forming) region is con¬ 
centrated at the center of the galaxy, and thus probes 
low values of v^. These options are consistent with the 
findings of state-of-the-art galaxy simulations which pre¬ 
dict t hat the gas in com pact SFGs reside in rotating 
disks (iZolotov et al.ll20T^ ). In the simulations, stars have 
°los U P 1-5 x larger depending on the projection of the 
gas disk. A clear prediction of these models is that the 
compact quiescent descendants should retain some rota¬ 
tion from its disky progenitors. 

GDN-8231 stresses the need for larger samples of com¬ 
pact SFGs with emission and absorption line kinemat¬ 
ics to quantify the effects predicted in the simulations. 
Those samples would allow us to study the dependence 
of ct los/°los w ith the viewing-angle and star-formation 
activity. Similarly, high spatial resolution imaging in the 
sub-millimeter with ALMA, or in the NIR with adaptive 
optics, can provide direct measurements of the size and 
location of the star-forming regions and, in some cases, 
resolved kinematics for the ionized gas. 
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TABLE 1 

Stellar and spectroscopic properties of the MOSFIRE sample 
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F160W 
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EW(Ha) 
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(3) 

(4) 

(5) 
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(9) 
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(13) 
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21.40 
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215±35 
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10.61 

0.87 

2.04 

295±110 

- 

aoi?j 

2 6" 1 ” 0 ' 7 

1.65±0.03 

1.31±0.27 


Note. — __ 

(а) Observed in lNewman et al.| j2010D ; IBelli et al.| d2014a|) , LRIS optical spectroscopy with a — 174km s — 1 . 

(I) General ID in the CANDELS iJ-band selected catalog in GOODS-N (Barro et al. in prep.) catalog. 

(2,3) R.A and Declination J2000. 

(4) Spectroscopic redshift. 

(5) H band (F160W) magnitude. 

(б) Stellar mass determined from SED fitting using [Bruzual Sz Chariot] J2003f) and a lChabrieif d2003T) IMF, see § |2.1| 

(7) Effective (half-light) radius (kpc) measured with GALFIT, see § 12.11 

(8) Sersic index measured with GALFIT, see § [2.1| 

(9) Integrated velocity dispersion (km s — 1 ) of the stars measured with PPXF dCappellari Sz Emsellem||2004l) , see §0 

(10) Integrated velocity dispersion (km s —1 ) of the gas measured from the emission line width (FWHM), see § 13.11 

(II) Equivalent Width of Ha (A) measured in the H band MOSFIRE spectra, see § |3.3| 

(12) Spectral index H5 a (A) measured in the Y band MOSFIRE spectra with PPXF, see § |3.3l 

(13) D4000 index measured from the SED including broad- and medium- band photometry as well as HST grism spectroscopy, see § |3.3| 

(14) Luminosity-weighted stellar age (Gyr) estimated from SED-fitting using the models of lPacifici et ah| d2012f) . 
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